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Microscopic Combinatorial approach is used to calculate the state and 
level densities with fixed exciton numbers, in some actinide nuclei. Deformed 
Saxon- Woods shell model was used as a basis from which all posible config- 
urations were generated. The pairing interaction was taken into account by 
applying the BCS theory to each configurations. Both the spin and parity 
distributions were obtained, considering the deformation of the nucleus. Rel- 
evance of the result to parity nonconservation studies involving epithermal 
neutrons on 238 U and 232 Th is discussed. 
PACS: 21.60w 

I. INTRODUCTION 

Adequate description of the nuclear level density of an excited nuclear state with a 
fixed number of quasiparticles is a basic ingredient of the statistical analysis of nuclear 
reactions. On the other hand, it is well known that the exciton dependent state density 
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used in the various pre-equilibrium models is lacking in considering the adequate particle- 
hole combination and its relation with the particular structure of the nuclei involved in the 
reactions. 

Analytical expressions for particle-hole excitation density can be obtained by means of 
methods of combinatorial or statistical mechanics in terms of the nointeracting particle 
model using the equidistant spacing approximation for the single-particle states. The main 
deficiency of the formulae commonly used to estimate exciton level densities comes from 
the assumption of equidistant spacings for the single particle levels U [@]. Usually only 
few-exciton configurations contribute significantly to pre-equilibrium emission and, in some 
cases , the population of low energy configurations is relevant as well. Therefore statistical 
approaches and equidistant spacings for the single particle levels do not seem to be adequate 
for pre-equilibrium calculations in some instances when shell structure effects are observed 
I- 

In a recent work [|]], two of the authors studied fundamental symmetry breaking on 
nuclear reaction and suggested that the observation of sign correlations in the longitudinal 
asymmetry of polarized neutron scattering from 232 Th at epithermal energies may constitute 
the first direct evidence for local 2p-lh doorways. The discovery of sign correlation in parity 
non- conserving (PNC) epithermal neutron-induced compound-nucleus reactions involving 
the heavy nucleus 232 Th have prompted intensive theoretical discussion concerning their 
origin. The great interest in the TRIPLE data 0, stems from the fact that the statistical 
theory of these reactions, although it predicts the possibility of large PNC at individual res- 
onances as the data also show, rules out any sign correlations in the longitudinal asymmetry, 
contrary to what the data show. More recently 0, data were obtained with smaller error 
bars for 232 Th and 238 U. Feshbach and collaborators started with the hypothesis that 
a very natural mechanism that could account for this sign correlations is to assume that 
the compound nuclear PNC process occurs through a single isolated p-wave local doorway 
which coupler to nearly s-resonances. In such studies, an important quantity which could 
determine how isolated the p-wave resonances are, is the 2p-lh density of states. It has been 
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a common practice to use the equidistant model to calculate the np-nh level densities, which 
is reasonable in spherical nuclei. No thorough investigation of the case of deformed nuclei 
is available. Sire nuclei such as U and Th are strongly deformed, one would be tempted to 
calculate the exciton level densities more accurately using a deformed mean field. 

Few years ago M. Herman and G. Reffb |7|, || evaluated a realistic few quasiparticle state 
densities with fixed exciton numbers using combinatorial calculations in spherical and some 
deformed nuclei. Their calculations were carried out in the basis of single particle orbitals 
derived from the harmonic oscillators defined by parameter due to Seeger and Howard [JTO] . 



Calculation of spin-dependence particle-hole level density were carried out only for spherical 
nuclei. For these nuclei the Gaussian- Wigner type spin factor was used. A direct test 
of the validity of this procedure is afforded by calculations made by Herman and Reffo. 
They found that the Gaussian formula works very well for configurations containing at 
least three excitons. However, in the case of deformed nuclei there are no comprehensive 
calculations with take into account the role of the collective rotational degrees of freedom 
for the evaluation of the particle-hole level density, with the question of the spin factor still 
open. 

It is generally assumed, that levels of different parities have an equiprobable condition. 
This is the Ericson's prediction || for a non-interacting many fermion system. The possibil- 
ity of considerable irregularities occurring in the distribution of positive-and negative-parity 
levels is demostrated in particular by the results of combinatorial calculations of level density 



carried out within the framework of the non-interacting particle model [^5[ as we H as those 
based on the more rigorous quasi-particle-phonon model, which takes into account the col- 
lective excitations of nuclei |26| |27]]. Previous combinatorial microscopical calculations |7j 
show that the results have evidently oscillatory character around the equal probability value. 
The fluctuation are observed throughout a wide energy range ( ~30 MeV),even though on 
the average the amplitude of the fluctuactions decreases with increasing excitation energy. 
The fluctuactions are reduced with the increase of the deformation, exciton number and 
mass of the nucleus. It is an indication that in most of the applications, the assumptions 



of equal parity distribution seems, to be justified. However, certain channels, particularly 
sensitive to the parity distribution (e.g. radiative channels), have to be treated with caution. 

Another interesting question is related to the deformation parameters used in the calcu- 
lation of single-particle spectra. Herman et al in || studied the effect of nuclear deformation 
on few-quasiparticle state-densities. It is shown that nuclear deformation tends to suppress 
strong fluctuations observed in the similar calculations performed in the space of spherical 
shell model basis vectors. The conclusion of Hermann and Reffo have the cualitative char- 
acter, since some fine details of particular structure of studied nuclei ( Al, 100 Mo, 170 Er, 
) are not taken in to account. In our view, the first step should be the right selection of 
deformation parameters with the calculations of single particle spectrum in the equilibrium 
deformation. The study of fine nuclear structure effects, of order of PNC need the most ac- 
curate calculation of particle hole level density, which is very sensitive to the single-particle 
spectrum. 

The paper is devoted to the development of a more accurate method for the calculation 
of particle-hole level densities for deformed nuclei by thoroughly considering the rotational 
degree of freedom. Particular effort is made in applying the method to actinide nuclei 
used in parity no-conservation studies induced by epithermal neutrons. The final intent is 
to evaluate the various uncertainties, constraints or difficulties related to the level density 
question and their impact on nuclear reaction predictions. 

In sect. 2 we outline very briefly the formalism with the description of the general method 
which improves our previous mode. The preliminary treatments. The problems of the spin 
and parity distributions is presented in sect. 3. In sect. 4 we discuss some results for two 
nuclei belonging to the actinide region. Finally, we present the main conclusions in sect. 5. 

II. PARTICLE-HOLE LEVEL DENSITY CALCULATION'S METHOD 

In order to calculate the 2p-lh doorway states densities coupled to total angular mo- 
mentum J we used the microscopic combinatorial method [|7| |§ JL1]] & t excitation energy E, 
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considering the axial deformation of the compound nucleus. The calculation was performed 
in the space of realistic shell-model single particle states, including pairing interaction in 
terms of the BCS theory. 

The first step to start the single-particle spectrum calculation is to fix the value of the 
equilibrium deformation of the nucleus using an optimal set of potential parameters JTJJ. In 



the shell model approach, based on mean-field potentials, the Strutinsky method is usually 
used fllBl . We have calculated the extreme points of the potential energy along the path 



from the equilibrium deformation to the more unstable form, using the BARRIER code 



| 12|1 which includes the Strutinsky prescription and the Pashkevich parametrization for the 



nuclear shape f2l |. The calculation was performed taking into account quadrupole and 



hexadecapole deformation parameters, which is a reasonable approximation to describe the 
equilibrium ground state. It is important to note that for more complex nuclear shapes, 



higher order deformation parameters should be included in the calculations ^0 . 

The other quantity which needs to be accurate by determined is the pairing strength 
parameter G njP . There are few methods commonly used to fix the value of G. Some authors 



take the pairing strength as a free parameter chosen to fit the calculated level 
densities to the experimental number of resonances at the neutron binding energy. This 
procedure certainly gives the appropriate level densities, but the pairing strength parameter 
derived this way is not related directly to the pairing interaction and may hide several 
inadequacies of the applied theory. We have chosen the method of relating the G parameter 
to the difference of masses of neighboring nuclei. 

In order to obtain intrinsic state densities all configurations with a specified number of 
particles and holes are generated within the assumed set of the single-particle levels by means 
of the permutation enumeration algorithm. For each of them a proper coupling of the spin 
projections is performed to obtain the nuclear states. The state density w(E,M) is found by 
counting states with the angular momentum projection M falling in the 0.5 MeV interval 
centered at the excitation energy. The pairing interaction was considered by applying the 
BCS model to each configuration so that a better understanding of the pairing correlations 



in few quasiparticle states could be achieved. 

We also used an improved formulae for the spin factor of particle-hole level density. 
Owing to the axial-symmetry of the ground states of these deformed nuclei we will use the 
spin factor formula which contains the rotational degree of freedom contribution, proposed by 
S. Bjornholm, A.Bohr and B.R.Mottelson in Ref The statistical microscopical method 



|TT| is used for to obtaining the parallel and perpendicular spin-cutoff factors in a consistent 



way. 



A. The single-particle spectra 

In the last few years, several calculations of single particle level schemes have been carried 
out both as a function of nucleon numbers and as a function of nuclear shape in connection 
with the studies of nuclear deformation potential energy surfaces. Simple global parameters 
for the underlying model potential have also been estimated on the basis of fits to nuclear 
ground state masses and fission barrier heights. Based on these and other single particle 
level schemes and the partition function approach, several numerical calculations of nuclear 
level densities have been carried out. 

In this work the nuclear shape parametrization is carried out using the Cassini 
parametrization [^TJ According to this the deforming shape (up to and beyond its sepa- 



ration into two fragments) can be conveniently described by the Cassini ovaloids figures 

Considering only axially symmetric nuclear shapes, the Cassini ovaloids are taken as the 
first approximation to the nuclear shape. The deviation from the ovaloid shape is given by 
an expansion into a series of Legendre polynomials. Geometrically, the family of Cassini 
ovaloids is defined by [?]: 

r 2 (z, e) = y/(a* + 4(cz)*) - (c 2 + z 2 - e 2 ). (1) 

In this equation, r and z are cylindrical coordinates; e is a dimensionless quantity such that 
c=eR 2 l ; c stands for the square distance from the focus of the Cassini ovaloids to the origin of 
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coordinates; and a is a dimensionless parameter which completely defines the shape taking 
into account volume conservation. 

In the plane containing the symmetry axis one can define a system of coordinates (R, x) 
such that the coordinate line R is constant. This is a Cassini ovaloid, where < R < oo 
and — 1 < x < 1. The (R,x) coordinates are related to the cylindrical ones (r, z) by the 
following equations 



R(z, r) = ^[{z 2 + r 2 ) 2 - 2eR% ■ {z 2 - r 2 ) + e 2 i$ (2) 

z 2 — r 2 — eR 



sign(z) 



1 + 



(3) 



R 2 (z,r) 

In this system of coordinates, the basic shape of the nucleus is described by these equa- 
tions, where R is constant, determining the Cassini ovaloids. Thus, the nuclear shape can 
be defined as a curve R(x) that does not intersect any straight line x = constant at more 
than one point. Accordingly, we can expand the function R(x) into multipoles, giving 

R(x) = Ro[l + (3 m Y m0 (x)}. (4) 

The set of parameters (e, (3) completely determine the nuclear shape. The details of this 
parametrization are given in ||12|| . As an example, we show in Figure 1 {e, a^} as functions of 
{^2,(^4}- As is clearly seen in this figure, it is difficult to establish an analytical connection 
between the two set of parameters. This relation was obtained by a least-square fit of the 
parameters /?2, ^4, for the harmonic spherical expansion, to our shapes described by the 
Cassini ovaloids. Using this figure it is possible to establish a connection between the two 
set of parameters to describe the same nuclear shape, but for more complex shapes more 
coefficients are needed in the harmonic spherical expansion. 

In order to obtain the single-particle energies and wave functions, the Hamiltonian ma- 
trix elements are calculated with the wave functions of a deformed axially symmetric os- 
cillator potential. The basis cut-off energy is determined in such a way that the negative 
energy eigenvalues of the Woods-Saxon potential do not change when adding more harmonic 
oscillator shells. 



As usual, the real potential V(r) is expected to follow roughly the density distribution. 
One of the most used radial dependences, is that of the Woods-Saxon potential, which 
takes into account the nuclear potential and density distribution. This potential involves 
the parameters Vo, r and a, describing the depth of the central potential, the radius and 
the diffuseness parameters, respectively. 

A definition of the radial dependence of the potential for a deformed nucleus, with an 



arbitrary shape of the surface, was proposed by Pashkevich ||21|| . According to Pashkevich, 
the nuclear potential is given by 

V(r,z,eJ) = (5) 

1 + exp 

where dist(r,z,e, p) is the distance between a point and the nuclear surface, a is the diffuseness 
parameter and Vo the depth of the potential well. 

The Woods-Saxon-type potentials, with the spin-orbit interaction proportional to the 
potential gradient, are the most appropriate from the physical point of view. The spin-orbit 
interaction is given by: 

V so (r, z, ej) = \ W(r, z, e, $) ■ (a x p) (6) 

where A denotes the strength of the spin-orbit potential and M is the nucleon mass. The 
vector-operator a stands for Pauli matrices and p is the linear momentum operator. 

The Coulomb potential is assumed to be that corresponding to the nuclear charge (Z — 
l)e, and uniformly distributed inside the nucleus.. In short, the depth of the central potential 
is parametrized as 

Oh 

V = V [1 + k'T z ], k' = ^ (7) 

with T z is the z-component of the isospin and k is an adjustable parameter. 

The single-particle spectra were obtained by means of the new version of the BAR- 
RIER code |12|| , based in a WSBETA code |T3|], using axial deformed Saxon- Woods nuclear 



potential well defined by parameters, due to Chepurnov |14 
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B. Microscopic Combinatorial approach for level densities with fixed exciton 

numbers 



The combinatorial method to obtain the level densities provides the possibility of direct 
counting of the levels with a fixed number of particle and holes. We start from the finite 
set of single particle states derived from the shell model with appropriately nuclear model 
potential, whose detailed description has been made above. Usually 100 neutron and 100 
proton orbitals were used in our calculation. It allows us to calculate the particle-hole state 
densities theory for deformed nuclei as Th 233 and U 238 , in energy range below 10 MeV. 

The pairing interaction was taken into account by applying the BCS theory to each 
configurations. All residual interactions but pairing are neglected. 

The total energy of each configuration is determined through the superconductivity the- 
ory. Configuration dependence is introduced into BCS theory by the blocking method as 
proposed by Wahlborn [0, to allow for blocking of more than are orbital. 

Accordingly for each generated configuration, a set of two BCS equations is solved 

n = 2 y, u 2 



- ■ e: [{si - \ s ? + a 2 ]- 1/2 



Ul = * 1- 77 2 J • 0) 



(8) 

G 

where 

1 [ 1 _ (f? ~ Xs 

2 V [(^-A s ) 2 + A 2 ] 7 

Here Si is the single particle energies, N stands for the number of paired nucleons, and A s 
and A s are the chemical potential and correlation function for a given configuration that 
are supposed to be determined. 

The total configuration energy, according to the BCS model 

E s = E * + 2E Ulei - ^ (10) 

The excitation energy is calculated in turn as the difference between the total energy of 
a configuration and the total energy of the ground state where the first summation included 
only blocked orbitals. 
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The particle-hole level density were obtained by means of the ICAR and CONV codes 
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C. Collective degrees of freedom. Spin dependence 



Collective phenomena in nuclei are receiving considerable attention in the analysis of 
spectroscopic data on the characteristics of low-lying levels. Various microscopic methods 



to consider the interrelationship between the collective excitations and the single-particle 
motion of nucleons in a self consistent nuclear potential. 

Strictly speaking, any separation of collective variables should be accompanied by a 
corresponding decrease in the number of internal degrees of freedom. But since collective 
motions are formed owing to deep-lying nucleons, while internal excitations are determined 
basically by the single-particle levels adjacent to the Fermi surface, exclusion of the extra 
degrees of freedom in the low-temperature region should not strongly affect the density of 
the internal excited states. Under these conditions, adiabatic consideration can be fully 
justified, at least as a first step in the analysis of the rotational increases, in the level density 
of nuclei. 

The impact of different models for spherical nuclei on the consideration of collective 
enhancement has been studied in Refs. pifl . Herman and Reffo J?]], have considered only 
internal degree of freedom at low energies without account of the rotation and the vibration 
enhancement can not be considered for these heavy nuclei The validity of the statistical law 
describing the spin distribution of nuclear levels must be reconsidered when applied to levels 
with fixed exciton numbers. 

For spherical nuclei the formula for the spin distribution function reads 



which is derived under the assumption of a Gaussian distribution of spin projections M. 
While it is very likely to be true when the number of levels is high enough, this assumption 



of describing the structure of collective levels |2?| |StJare also widely used at present 




2J + 1 



e -[(J+i) 2 /2 CT 2 ] 



(11) 
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may not hold for levels with low exciton numbers, for which the density of states is too low 
for statistical treatment, to be justified. 

In the case of spherical nuclei Herman and Reffo [|7| found that Eq. (5) is valid for levels 
containing at least four and, to some extent, also three excitons. 

The collective contribution to the level density of a deformed nucleus is defined by the 
symmetry order of nuclear deformation. For deformed nuclei the spectrum of its energy 
states will be determined not only by the internal excitations but also by the rotation of 
the nucleus as a whole. This rotation may lead to a considerable increase in the density of 
nuclear levels. Since the deformation of the nuclear potential removes the degeneracy of the 
basis vectors belonging to the same spin multiplet, we are no longer able to obtain the spin 
distribution of nuclear levels, according to Eq. (5). 

Reffo et al || have studied the influence of nuclear deformation on the distributions of 
quasiparticle states. It is shown that deformation tends to suppress the strong fluctuations 
observed in the similar calculations performed in the space of spherical shell-model basis 
vectors. In a deformed nucleus, each intrinsic state gives rise to a rotational band on the 
total level spectrum, for a given angular momentum. It is therefore obtained by summing 
over a set of decomposition of the level spectrum, as for a spherical system. The experimental 
level data of the nuclei under study show that this prediction is well established. In that 
case, the following expression has been used for level-density calculation of axially-symmetric 
nuclei 



-Jji+2) 
— 



(12) 



R(J) = *±L. e 

v ' Vo7r0-|| 

P2 P ih (U, J) = uj 2p ih(U) . R(J) ■ 

In this formula the contribution of rotational states is taken into account. Here uj^pih (U) is 
the microscopic level density of 2p-lh states. 



The spin cut off parameters aj_ and <7m are calculated in the following way |L9 



ai = ^ af { =Q- 2 gt- (13) 
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where t is the nuclear temperature and g is the single-particle level density near the Fermi 
energy, and fl~ 2 is the value of the average single-particle square projection on the symmetry 
axis of deformed nucleus. 

The energy dependence of the moment of inertia 3j_ is approximated in the following 

way: 



= { Uc " tJ (14) 

cv jj \ tj 

~^>rtg u ^ ^ crit 

So is the moment of inertia in the ground state, U cr # is the maximum value of the transition 
energy from the superfluid to the normal state for neutron and proton system, and Qyj g is 
the rigid body moment of inertia of the nucleus (^s r ig = ^Mr^A^). 



D. The parity distribution 

According to Ericson |§, even and odd parity levels contribute equally to the level density. 
Being aware of the fact that the equal parity distribution is very questionable for levels with 
a fixed number of exciton, we have performed combinatorial calculations to investigate this 
problem. This distribution at a fixed value of energy can be described in terms of the 
asymmetry ratio: 

A{U) = n+(U) + n-(U) (15) 

where N 71 " are the number of levels with parity n — ±1. 



III. CALCULATIONS AND RESULTS 



Based on the above described procedure, we consider in the following 238 U and 232 Th. 
Both nuclei are well deformed in their ground states and therefore, when they interact with 
epithermal neutrons the compound nucleus acquire reasonably high deformations. For this 
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reason, these two nuclei constitute a natural test of the procedure proposed in the present 
work. 

Figures la and lb show the potential energy surface calculated in the frame of the 



Strutinsky's Method |I5[] and the corresponding ground state deformation points for each 
case which are in agreement with the experimental data. As already mentioned, a com- 
prehensive description of the level densities should include the characteristic of the nuclear 
configuration for different deformation parameters. In this sense, we show also other extreme 
points which define a possible path in the multidimensional space of the deformation de- 
grees of freedom. This path could describe a decay-process of the compound system from its 
ground state and therefore has to be considered in calculations of cross sections in particular 
in reactions induced by epithermal neutrons. 

Usually in calculations dealing with such highly deformed and excited compound nuclei 
where a drastic change of the nuclear shape ensues, one relies on extrapolations from the 
ground state structure properties. Of course, this represents a very crude procedure because 
in these regions the shell closure conditions are changing and one has to take into account 
the role of higher multipole coefficients to describe the nuclear shape. 

The Chepurnov parameters, and the obtained deformation parameters allow us to repro- 
duce the experimental low-lying quasi particle states of the nuclei and the ground state of 
neighboring nuclei. 

The extremal points in Potential Energy Surface of Th 233 are given in Table 1. 
The single particle level density is a basic ingredient in the calculations, because the 
single particle spectra depend on the symmetry properties of the potential well. In nuclei 
deformed in their ground states the influence of the second minimum on the level density will 
be sharply decreased because of a smaller difference between deformation in two minima. 
This situation is showed in Table 2. 

The calculations of the particle-hole level densities were carried out using the modified 
version of ICAR code ]T7| . The results of the calculation for 2p-lh states are shown in 



Fig. 2 for positive parity in comparison with the Williams formula. The solid line represents 
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the Williams formula, while the dashed lines corresponds to calculation obtained with our 
method using the Cherpunov parameters of Wood-Saxon nuclear potential. The agreement 
is reasonable, though the Williams formula represents an average of the combinatorial level 
density. 

The pairing strength parameter G is calculated from the difference of masses of neigh- 
boring nuclei as proposed by M. Herman @. The values obtained for 232 Th and 238 U are 
given in Table 3. 

In Fig. 3 we show the comparison of the behavior of the spin cutoff parameters o\ and 
a 2 as function of energy, for Th 233 and U 239 . The solid lines represent the results for the 
Th isotope using the approach described in eq.(13), while the dashed lines correspond to 
the calculations for the U isotope. In the case of a± the figure shows a similar behavior 
for both isotopes, which reflects the similarity of the moments of inertia calculated in our 
approach; this results from almost equal shapes found in the ground state. In the case of cry , 
the calculations exhibit marked differences for the two isotopes. This fact reveals the need 
of taking into account the single-particle spectrum for each particular isotope, since for crjj 
this plays a crucial role. 

In Fig. 4 we show the behaviour of the spin factor distribution as a function of energy, for 
Th 233 in the spherical and the deformed approaches. Note that the statistical assumption of 
a Gaussian distribution with spherical symmetry is not always justified, because for the cases 
of low J (in Fig. 4 we show for J=|, see also Figs. 5 and 6) it underestimates the calculations, 
which has an important influence on analyses such as those performed, in ||. In Fig.5 we 
compare the influence of the spin cutoff parameters on the surface of the distribution R(E,J), 
in comparison with the spherical formulae using eq.(ll) for low energy and J- values. The 
absolute values differ by nearly two orders of magnitude for the spherical and deformed 
cases. 

For small p-h excitations the exact combinatorial solutions described in Sect. II should 
be used. To take into account the rotational degrees of freedom, the rotational bands are 
built on each quasiparticle state, using the average values of the rotational constants A r and 

14 



E(J, K*) = E(K W ) + [j ( J + 1) - K 2 ] A r + B r ^J+ (-1) 8 K i 

The low energy region is interesting from the point of view of its difference from the 
region where statistics is applicable. The number of levels is small and their energy and 
quantum number distributions may have strong fluctuations. 

For large p-h excitations, the p-h level densities with spin distribution should be used with 
the approach described in Eq.(12). The width of the spin distribution can be determined 
using Eqs.(13) and (14). This procedure provides a convenient way for the analysis of 
deformed nuclei. 

The results of the microscopic calculations of two-particle-one hole level densities for J=|, 
in Th 233 and U 238 , are shown in Fig. 6. The impact of collective rotational enhancement 
in the level density can result values 3 times large than those from the Williams formula 
using the standard parameters for these nuclei. This procedure can be useful in order to 
reevaluate the parameters of the Williams formula. 

Finally, in Fig. 7 the trend of the fraction of positive parity levels for the analyzed nuclei 
as a function of excitation energy (up to 15 MeV) is plotted. Here, this fraction is given by 
the parity asymmetry defined by Eq.(15) as a function of U, where N ± (f/) is calculated by a 
direct counting of the levels. The regions of lower energies show the greatest fluctuations in A 
(eq. 15); although these fluctuations persist at higher energies, a slow smooth out is observed. 
This observed strong oscillatory character, not predicted by the Ericson hypothesis, must 
be taken into account in studies of PNC, for which the fine structure effects are important. 



IV. CONCLUSION AND FINAL REMARKS 

We have presented a realistic few quasiparticle level density approach for deformed nuclei, 
which should be preferred to the more phenomenological expression commonly used. The 
more realistic inclusion of the rotational degree of freedom improves the calculations in 
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the region of the actinide nuclei, where deformation is not negligible even in the ground 
state. The analysis of nuclear reactions in the region of the actinide nuclei, with epithermal 
neutrons, requires the consideration of nuclear deformations where the amplification of the 
collective movement must be taken into account. The deformation introduces significant 
fluctuations in the level densities. This prevents the use of a single analytical expression 
able to reproduce these densities in different mass regions. In this regard, a judicious analysis 
of the extreme points becomes essential. When one needs to analyze the shape of the nucleus 
corresponding to the external saddle point, it is necessary to use a procedure like the one here 
proposed, because the nuclear deformation for a fixed configuration of particles and holes 
neither suppresses the fluctuations that takes place in the level densities, nor the parity 
distribution. The importance of our calculation for an eventual quantitative understanding 
of the mechanism of PNC in epithermal neutron induced compound reactions using actinide 
nuclei is discussed. 

V. ACKNOWLEDGMENTS 

This work was partially supported by the CNPq-CLAF and FAPESP, and in part by 
funds provided by the U.S. Department of energy (D.O.E.) under Cooperative Agreement 
# DE-FC02-94ER40818. 



16 



VI. REFERENCES 



[1] E.Gadioli and P.E.Hodgson, in Pre-Equilibrium Nuclear Reactions, Oxford Studies in Nuclear 
Physics Vol.15, Oxford Science Publications (1992). 

[2] W.Scobel, M.Blann, T.T.Komoto, M.Trabandt, S.M.Grimes, L.F.Hansen, C.Wong and 
B.A.Pohl, Phys.Rev.C30, 1480 (1984) 

[3] S.Shlomo, Ye.A.Bogila, V.M.Kolomietz and A.I.Sanzhur, Z.Phys.A353, 27 (1995). 

[4] H. Feshbach,M.S. Hussein and A.K.Kerman, Fundamental symmetry breaking studies with 
nuclear reactions, in International Workshop on Parity and Time reversal Violation in Com- 
pound Nuclear States and Related Topics, edited by N. Auerbach and J.D. Bowman, World 
Scientific, Singapure, P. 157 (1996). 

[5] C.M.Frankle et al ,Phys.Rev Lett 67,564 (1991) 

[6] F. Mitchel, J.D. Bowman and H. A. Weidenmiiller, Rev. Mod. Phys. 71, 445 (1999) 

[7] M. Herman and G. Reffo, Phys. Rev. C 36 N4, 1546 (1987). 

[8] M. Herman, G. Reffo and R.A. Rego, Phys. Rev. C 37 N2, 797 (1988). 

[9] Ericson T., Ann. of Phys. 23 (1963) 390. 

[10] PA. Seeger and W.M. Howard, Nucl. Phys. A238, 491(1975) 

[11] F. Garcia et al., J. Phys. G19, 2165 (1993). 

[12] F. Garcia, O. Rodriguez, J.D.T. Arruda-Neto,E. Garrote, J. Mesa,F. Guzman and R. Capote, 
BARRIER code : accepted for publication in Comp. Phys. Comm. (1999). 

[13] S. Cwiok et al., Com. Phys. Comm. 46, 379 (1987). 



17 



[14] V.A. Chepurnov, Sov.J.Nucl.Phys,6, 696 (1968). 

[15] V. Strutinsky, Nucl. Phys. A95, 420 (1967). 

[16] Wahlborn, Nucl. Phys. 37, 554 (1962). 

[17] M. Herman and G. Reffo, Com. Phys. Comm.47,103 (1987) 

[18] P. Moller, J.R.Nix, W.D.Myers and W.S.Swiatecki, Atomic Data and Nuclear Data Tables 
59,185(1995) 

[19] F.Garcia, O. Rodriguez, V.A.Rubchenya, E.Garrote, Comp. Phys. Comm. 86 (1995) 129. 
[20] J.D.T.Arruda-Neto et al, Phys.Rev.C55,(1997). 
[21] V.V.Pashkevich, Nucl.Phys.A169, 275 (1971). 

[22] G.Maino, E.Menapace and A.Ventura, Nuovo Cimento 57A, 427 (1980). 
[23] V.Benzi, G.Maino and E.Menapace, Nuovo Cimento 66A, 1 (1981). 

[24] S. Bjornholm, A. Bohr and B.R. Mottelson, ,Proceedings of the Third IAEA Symposium on 
the Physics and Chemistry of Fission, IAEA in Rochester New York, Vol I, 367 (1973) 

[25] M.Hillman and J.Grover, Phys.Rev.,Vol.l85, 1303 (1969). 

[26] V.G.Soloviev, Ch.Stoyanov and A.I.Vdovin, Nucl.Phys.A 224, 411 (1974). 

[27] A.I.Vdovin et al, Sov.J.Part and Nucl.,Vol.7, N4, 952 (1976). 

[28] A.B.Migdal, The theory of finite Fermi systems and the properties of atomic nuclei, Nauka 
Press,Moscow(1965) . 

[29] V.G.Soloviev, Theory of complex nuclei, Ed. Mir Moscow (1971) 

[30] V.G.Soloviev, Theory of Atomic nuclei: Quasiparticles and phonons, Ed. Mir Moscow (1989) 

[31] J.M.Eisenberg and W.Greiner, Nuclear Models, Vol.1: Collective and Single-Particle Phe- 
nomena, Elsevier Science Publishers B.V.,New York (1987) 

18 



[32] A.Borh and B.R. Mottelson, Nuclear Structure, Vol. 1,2 W.A. Benjamin Inc. N.Y.-Amsterdam 
(1969) 



19 



VII. TABLE CAPTIONS 



Table 1. Deformation parameters of extremal points in potential energy surface of 
Table 2. Single-particle level densities of extremal points in 233 Th. 
Table 3. Pairing constants of 233 Th and 238 U. 
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VIII. FIGURE CAPTIONS 



Figure la. 233 Th Surface Potential Energy. See text for details. 
Figure lb. 239 U Surface Potential Energy. See text for details. 

Figure 2. Microscopic combinatorial 2p-lh level density corresponding to positive parity 
in comparison with the William formula. See text for details. 

Figure 3. Spin cuttof parameters vs. excitation energy. See text for details. 

Figure 4. The level density spin distribution factor R(J) for 233 Th vs. excitation energy 

Figure 5. Two-dimensional Plot of R(J) for 233 Th. See text for details. 

Figure 6. The 2p-lh level density for 233 Th and 239 U. The Microscopic Combinatorial 
Method results are compared with the equidistant Williams formula. 

Figure 7. Percentage of positive parity levels in 233 Th and 239 U. See text for details. 
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IX. TABLES 



Table 1. 



Def. param 


1st Min. 


1 Saddle 


2d min. 


3th min. 


2 Saddle 


3 Saddle 


e 


0.22 


0.383 


0.49 


0.7213 


0.7044 


0.755 


a 3 


0.0 


0.0 


0.0 


±0.1179 


0.0 


±0.087 


0:4 


0.071 


-0.0827 


0.0275 


0.0386 


0.0386 


0.0386 



Table 2. 



Ext. point 









G.st 


3.651 


2.693 


6.345 


2 Min 


3.705 


2.669 


6.374 


2d Sadd Asymm 


3.77 


2.58 


6.35 


2d Sadd Symm 


4.034 


3.069 


7.103 


3rd Min 


3.661 


2.609 


6.27 



Table 3. 





Chepurnov parameters 


Optimal parameters 


Th 233 


20.25 


18.5 


TJ239 


19.9 


18.0 
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Fig 1a. 
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Fig.3 




Fig. 4. 



1.00 



J=1/2 



Spherical 

— — Deformed 



0.10 



0.01 



0.00 



0.00 



4.00 8.00 12.00 16.00 

Excitation Energy (MeV) 



20.C 



Fig5. 






0.00 2.00 4.00 6.00 8.00 

Excitation energy (MeV) 



Fig. 6. 



V) 

a> 
> 
a> 



a> 
> 
a> 



re 

Q. 

a> 
> 



0.52 



Uranium 239 



■5 0.48 
o 

Q_ 




0.00 4.00 8.00 12.00 16.00 

Excitation energy (MeV) 



V) 

a> 
> 
a> 



re 



v> 

a> 
> 
a> 



re 

Q. 

a> 
> 

o 

Q. 



0.52 



Thorium 233 



0.48 




0.00 4.00 8.00 12.00 16.00 

Excitation energy (MeV) 



Fig.7 



